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ABSTRACT 


A discussion of the generation of the Best Estimate Trajectory (BET) 
of the first NASA Space Shuttle Orbiter entry flight (STS-1) as reported by 
Compton, et al. , in Reference 1 is presented. This work was sponsored 
by NASA LaRC under Contract No. KAS1-16087 to the Analytical Mechanics 

Associates, Inc. The BET defines a time history of the state, attitude, 

and (combined with the best available atmosphere as defined by the Langley 
Atmosphere Information Retrieval System (LAIRS) ) atmospheric relative 
parameters throughout the Shuttle entry from an altitude of approximately 
183 km to rollout on Runway 23 on the Roger's dry lake bed at Edwards Air 
Force Base. The inertial parameters were estimated utilizing a weighted 
least squares batch filter algorithm. Spacecraft angular rate and acceler- 
ation data derived from the Inertial Measurement Unit (IMU) were utilized 
to predict the state and attitude which was constrained in a weighted least 
squares process to fit external tracking data consisting of ground based 
S-band and C-band data. In addition, refined spacecraft altitude and velocity 
during and post rollout were obtained by processing artificial altimeter and 
Doppler data. 

Appendix A is presented to provide for a general discussion of the BET 
generation process. This' includes both software and data interface discussions 
as well as a definition of the variables and coordinate systems utilized. STS-1 
mission peculiar inputs are summarized in Appendix ®, though the report 
contains tables and figures which show the more relevant results, it is virtually 
impossible to present all the information in this form! Thus, Appendix C is 
included which provides a listing of the contents of the actual BET. 



I. Introduction 


The completion of the first successful flight of the Space Shuttle Columbia 
on April 14, 1981 opened a new era in NASA's manned spaceflight. Researchers 
at the NASA Langley Research Center, as well as others throughout the aero- 
space community, have proposed use of the Shuttle as a research vehicle for 
postflight aerodynamic and aerothermodynamic investigations (References 2, 

, and 4). The best postflight trajectory and atmospheric information is a 
necessary input for such investigations as the Aerodynamic Coefficient Mea- 
surement Experiment (ACME). Development of the best available atmosphere 
based on models as well as meteorological measurements is discussed in 
Reference 5. This report discusses the generation of the required trajectory 
information using the methods discussed by Compton, et al (Ref. 1,6). The 
process is functionally presented as Appendix A of this report in terms of a 

software overview and the required pre-processing of both the observational 
and dynamic data. 

AMA, Inc. , under NAS1-16087, is responsible for this postflight tra- 
jectory reconstruction, as well as generation of the final product for use by 
the user community. The reconstructed trajectory, based on onboard mea- 
surements of the spacecraft dynamics and ground based radar tracking, is 
necessarily an inertial product. To satisfy the total requirements of the 
aerodynamic and aerothermodynamic researchers, the final product (Ref. 7) 
merges the inertial reconstructed entry history with the best available atmos- 
pheric data. This product includes computation of the important atmospheric 
relative parameters as well as first order estimates of the flight derived total 
aerodynamic coefficients. 

Section II presents a procedural discussion and includes an overview of 
the tracking coverages for STS-1. Mission specific input data are presented 
as Appendix B. Results are presented in Section in. Section IV summarizes 
these results and presents conclusions. Finally, a listing of the STS-1 BET 
parameters is presented as Appendix C. 
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II. Procedural Discussion 
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n. 1 Mission and Spacecraft Specific Data 

There are numerous flight-dependent inputs required by the 
various elements of the entry reconstruction software, ENTREE (Ref. 8). 

These are given in Appendix B. Tracking station locations, acronyms, and 
refraction constants are given in Table B-l. These data were obtained from 
the mission software data base, Revision G.02 (Ref. 10). The required IMU 
attitude transformation matrices are given in Table B-2. These data were 
obtained from the Johnson Space Center and Ref. 9. Assumed a priori 
parameter uncertainties are given in Table B-3. Planet model parameters, 
Runway 23 locations, IMU locations with respect to the Shuttle center-of- 
gravity and Shuttle mass properties and aerodynamic reference values are 
presented in Table B-4. 

II. 2 Initial Condition State Vector 

Initial position and velocity estimates in Cartesian Mean of 
1950 (M50) coordinates were provided by the Math Physics Branch at JSC. 

This state vector was the real-time Guam tracking pass solution and was valid 
at 17 h 42™30 S GMT on April 14, 1981. Since the time was very close to Guam 
Acquisition of Signal (AOS), it was chosen as the epoch (63750? 0 from midnight, 
day of entry) for the STS-1 BET. The 6 -element state was transformed to 
ENTREE input coordinates (spherical, Earth-fixed, Earth true equator of 
date) using standard formulas. Figures A-3a and A-3b in Appendix A define 
the ENTREE variables of interest. Initial attitude estimates (one per each 
IMU) were obtained using the attitude transformation matrices given in 
Appendix B, the 6-element state, and the interpolated platform to outer roll 
quaternions (at the state vector epoch) from the telemetry tape. The resulting 
start vector conditions are shown in Table n-1. Note the consistency in attitude 
estimates among the IMUs. 

II. 3 Dynamic Data 

Dynamic data, which consists of measured spacecraft angular 
rates and linear accelerations, are required for the BET generation. This 
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requirement was satisfied by the IMU measurements. A performance eval- 
uation among the three onboard IMUs (Ref. 11) showed very good consistency 
in their respective measurements. Based on this analysis and other com- 
parisons of the IMU derived dynamic data, no "preferred" IMU could be 
determined. Since IMU2 had shown perhaps the best trajectory prediction 
capability (using initial condition estimates obtained from JSC), it was 
selected as the primary dynamic data source for BET development. How- 
ever, as will be shown in Section III, very good trajectory solutions were also 
obtained using IMU1 and IMU3. 

Essentially continuous measurements, i.e. , no major data gaps, were 
obtained from each of the IMUs. IMU data covering the entire entry from the 
Guam AOS to approximately 17 S after vehicle stop were used. The only cor- 
rection made to the "raw" data was a 0. 007 sec adjustment to account for the 

spacecraft clock lagging the station clocks. This clock offset was provided 
by the JSC. 

Figures Il-la through II-lc show the dynamics experienced by the space- 
craft during the STS-1 entry flight. Plotted are the body axis components of 
the angular rates (Fig. n-la), the linear accelerations (Fig. II- lb) and the 
angular accelerations (Fig. II-lc). These data were derived from the 1 Hz 
(nominally) IMU2 measurements using the methods described in Appendix A. 
The spacecraft rates and accelerations in the platform frame were rotated to 
the body axes and translated to the vehicle center-of-gravity. Angular accel- 
erations were obtained by numerically differentiating the angular rate data. 

II. 4 Tracking Data 

Radar tracking data from the Guam S-band station and eight (8) 
California C-band stations were used in reconstructing the STS-1 entry tra- 
jectory. Appendix B contains a list of the station acronyms, locations, and 
refraction constants. Appendix A describes the pre-processing required, in 
general, pre-processing was very straightforward and consisted primarily 
of reordering and units conversions. However, the Guam high speed S-band 
data obtained from GSFC required time-tag corrections. According to GSFC, 
this problem is unique to playback data and can be expected on subsequent 
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flights. The time-tag corrections were made using low speed real time 
listings obtained from both GSFC and JSC. The adjustments made are 
given below in terms of GMT time on April 14, 1981 and also, in paren- 
theses, the time from the BET reference epoch. 

• Range, Doppler from 17:44:16. 3 (106?3) to the end of the pass 

s 

were time-shifted earlier by 0. 1 

• X, Y-angles from 17:42:18 (-12?0) to 17:44:16.3 (106?3) were 

Q S 

time-shifted earlier by 0.1 and from 17:44:16.3 (106.3) to the 

s 

end of the pass were time-shifted earlier by 0.2 

Fig. II-2 presents the complete STS-1 entry ground track (~ 40 min) 
overlaid on a geographical map segment. Also indicated are the tracking 
sites and approximate spacecraft altitudes at 500 sec increments along the 
track. 

Tables TI-2 and II-3 together with Figs. II-3a through n-3c illustrate 
the detailed tracking coverage. Table n-2 is a sequence of events for the 
trackers and shows acquisition of signal (AOS), loss of signal (LOS), and 
maximum elevation during the pass. Also, approximate observations are 
given at the specific times for information. In the case for the S-band station 
(GWMS), derived elevation data are shown. Table II-3 indicates the actual 
data arc processed for each tracker, subject to the processing constraints 
(principally elevation angle cutoff) used. 

Figure n-3 presents the station coverage during each of the three main 
entry segments. The coverage for each station is shown by "rays" from the 
station to the ground track. Coverages indicated are the actual arcs pro- 
cessed (Table n-3). Also, for better illustration, only one station from the 
Vandenberg and Pt. Pillar complexes are shown. Coverage for the other 
stations in these complexes is similar. 

The limited upper altitude coverage and the importance of the Guam pass 
are shown in Fig. n-3a. In time and altitude, the Guam pass covers approxi- 
mately three (3) minutes and an altitude range from — 183 km to ~ 145 km. 

The C-band stations were not acquired until approximately 21 minutes after 
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Guam LOS at an altitude of ~ 55 km. (The first C-band measurement pro- 
cessed was at 1577. 0 corresponding to an altitude of ~ 50 km). Fig. H-3b 
indicated considerable overlapping C-band coverage for approximately six(6) 
minutes over the altitude range from ~ 50 km to ~ 23 km. Fig. n-3c shows 
that during the last 6 minutes of the entry, from h~ 23 km to h~ . 06 km, 
only Edwards and Dryden coverage was available. Dryden tracking lasted 

until main gear touchdown, whereas Edwards coverage ended about 17. S 0 
earlier. 

In summary, for a 40 minute entry, radar tracking data processed were: 
(1) approximately three(3) minutes of high altitude coverage (183 km to 145 km) 
from Guam; (2) approximately six(6) minutes of 8-station overlapping C-band 
coverage (50 km to 23 km); (3) approximately five(5) minutes of the dual station 
coverage from approach to landing (23 km to . 06 km). 

All tracking data were processed at a 2 second data rate. A five(5) 
degree elevation angle cutoff constraint was used. An exception to this was 
the Dryden and Edwards Range and Azimuth data to enable better coverage at 
touchdown. The assumed data accuracies were based on preflight specifications 
and the actual scatter in fit residuals during processing. Assumed S-band 
accuracies were 1.5 m for Range; 0.3 Hz (~ 20 mm/sec) for Doppler; 0.2 mrad 
for both X and Y-angles. Those for C-band were: 9m for Range; 0. 2 mrad 
for both Azimuth and Elevation angles. S-band X-angles were not processed 
when Y-angle measurements exceeded 70 degrees because of known X-angle 
inaccuracies in this region. In addition, C-band angles were not processed 
when the spacecraft was near zenith over Edwards and Dryden. All radar 
measurements, except C-band Azimuth, were corrected for atmospheric re- 
fraction using the algorithm given in Ref. 12. The modulus of refraction at 
each station was the mean monthly value for April as shown on Table B-2. 
Atmospheric scale heights were obtained using the algorithms of Ref. 12. 
Tracking observations were also corrected for the light-time delay using 
extensions of the procedures described in Ref. 13. 
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IT. 5 Other Observations 


In addition to the C-band and S-band tracking data, two types 
of pseudo data were processed during and post rollout on the dry lake bed. 
During rollout, the vehicle c. g. is known to be about 4. 8768 m above ground 
level, within ± 1 m due to strut deflections resulting from various aerody- 
namic and wheel brake loads acting on the vehicle. Thus, pseudo altimeter 
observations of 4. 8768 m were processed every second from t = 2318. 0 
(following nosewheel touchdown) through the end of the estimation run at 
t = 2384? 0 (16 seconds following vehicle stop). The altimeter data were 
weighted to an assumed 1 m (Iff ) accuracy. In addition, beginning at 
t = 2370?0, pseudo Doppler data consisting of 0.0 Hz (hull) observations 
were processed 1 per second from 3 ficticious S-band stations located 609. 6 m 
to the North, East, and below the vehicle stop position. The pseudo Doppler 
data were weighted to an assumed accuracy of 0. 1 Hz (Iff ). Inclusion of 
these pseudo measurements, which were based on known terminal flight con- 
ditions, rectified the BET trajectory to eliminate approximate errors of 

0. 4 mps and 17 m velocity and altitude, respectively, during and post rollout. 

n. 6 Solution Parameter Selection 

During the reconstruction process, in addition to solving for 
the required spacecraft position, velocity and attitude, inclusion of both 
dynamic and observational parameters as solution parameters in the estimation 
was considered. Although many sets of these "extended solve-for parameters" 
were studied, the final BET included only six: 3 IMU gyro drifts, and 3 IMlT 
accelerometer scale factors. Ideally, if the dynamic and observational instru- 
ments were perfect, the BET could be determined via a state-only solution, 

1. e. , position, velocity and attitude at epoch. However, the total weighted 
root mean square (RMSW) of the tracking residuals for a state-only solution 
was 2. 2. In other words, the overall fit was 2.2 times the assumed 1 a 
accuracy of the tracking measurements. Although the state only solution 
provided reasonable initial and terminal state vectors, additional parameters 
were included in the solution set to improve the fit to the tracking data and 
obtain a better entry trajectory. 
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Many factors influenced the final state vector size selection. First, it 
was believed that solving for observation related biases would not really im- 
prove the estimation accuracy though the data fit might appear to be better in 
the sense that the mean errors were reduced. It v/as felt that the best way to 
account for any potential measurement related error source was to process 
the data from all^ available stations, thus, in effect, averaging the errors, if 
any. Thus, the final BET was determined from the uncorrected tracking 
data. 

Pre-mission simulations had shown that (1) center of gravity position 
errors many times larger than the Uncertainty associated with the advertised 
c. g. location had a very small effect on the ensuing estimation accuracy, and 
(2) with the tracking data accuracies available, little if any c.g. location 
information could be extracted from the data arcs. Hence, center -of-gravity 
errors were not solved for. 

Early studies were done with various combinations of eighteen(18) poten- 
tial IMU error sources in ENTREE: accelerometer biases (3), accelerometer 
scale factors (3), gyro drift biases (3) , and g-sensitive gyro drift biases (9). 
Note that since only body to actual platform attitude information is necessary 
to derive the dynamic data for ENTREE, any initial IMU misalignments re- 
sulting from the pre-deorbit star tracker alignment need not be modeled or 
solved for. 

With the previously mentioned 18 instrument parameters included in the 
solution set, the RMSW was reduced to 1. 02. However, removing the 9 
g-sensitive terms hardly degraded the fit, i. e. , the RMSW increased to 1. 05. 
Also, the dependence on a priori was reduced when g-sensitive terms were 
eliminated. Furthermore, based on conversations with JSC flight controllers 
who indicated that a successful pre-deorbit accelerometer calibration had 
transpired, and based on IMU comparisons (ref. 11) which indicated acceler- 
ometer bias errors on the order of only 10 yg, the 3 accelerometer bias 
parameters were also removed from the solution set. This left the 3 acceler- 
ometer scale factor errors, and the 3 gyro drift bias errors in the extended 
solution set of the final BET. 
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Cartesian M50 


ENTREE Coordinates 


-2370-97465 km 
-6113.30502 km 
+ 226.76197 km 
+5.645572676 km/sec 
-1.843230515 km/sec 
+5. 008001519 km/ sec 


7.4108907 km/sec 


-1. 1568500 deg 


47.213181 deg 


182. 76046 km 


1.9311855 deg 
140. 76250 deg 
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* 

Time 

Site 

Event 

Range 

Azimuth 

Elevation 

X-Angle 

Y-Angle 




(km) 

(deg) 

/Hpcr\ 


-..(dee) _J 






.Uis&J. . .. 

— iS*£KX 

0 

GWMS 

AOS 

1341 


1. 7 + 

-83.6 

-70.5 

155 

GWMS 

max elevation 

671 


11.2 + 

71.9 

-51.0 

313 

GWMS 

LOS 

1280 


0.3 + 

88.1 

12.9 

1522 

VDBC 

AOS 

579 

284.2 

2.7 

> 

r/A 

1534 

VDFC 

AOS 

549 

286.0 

3.1 



1535 

SNIC 

AOS 

701 

296.1 

1.2 



1574 

FRCC 

AOS 

640 

280.9 

1.7 



1577 

VDSC 

AOS 

427 

291.4 

5.0 



1583 

EAFC 

AOS 

810 

281.4 

1.9 



1632 

PPTC 

AOS 

177 

200.1 

14.9 



1650 

PPTC 

max elevation 

165 

184.7 

15.6 



1715 

PTPC 

AOS, max elevation 

216 

142.9 

9.9 



1767 

VDBC 

max elevation 

125 

20.1 

16.5 


| 

1768 

VDFC 

max elevation 

131 

19.6 

15.2 



1769 

VDSC 

max elevation 

131 

19.9 

15.2 



1834 

SNIC 

max elevation 

223 

5.6 

7.5 



1893 

PPTC 

LOS 

427 

127.1 

1.8 



1910 

PTPC 

LOS 

457 

126.1 

1.9 



2012 

EAFC 

max elevation 

17 

18.3 

84.1 



2018 

FRCC 

max elevation 

16 

18.3 

82.8 



2137 

VDBC 

LOS 

274 

82.8 

1.0 



2149 

VDFC 

LOS 

274 

82.5 

-2.2 



2156 

VDSC 

LOS 

274 

77.5 

-1.1 



2162 

SNIC 

LOS 

262 

41.9 

1.8 



2281 

EAFC 

LOS 

12 

87.8 

-0.6 



2305 

FRCC 

LOS 

7 

90.7 

-1.0 

1 ' 

J 


* J 1-1 

Time in seconds from epoch 104 a i7 42 m 30 S 
Derived for information only 

TABLE II- 2 

STS-1 C-band and S-band Sequence of Events 
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deg /sec 



Figure n-la. STS-l body axis angular rate history derived from IMU2 measurements 




Figure n-lb. STS-1 body axis acceleration history derived from IMU2 measurements 





angular acceleration history derived from IMU2 measurements 







TIME TICS* 0. TO 0000. AT 



Figure n-2. STS-1 entry ground track 










LATITUDE (DEG) 





(c) Final Approach and Landing 
Figure II-3. Detailed tracking coverage geometry for STS-1 
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III. Results 


Though most of the results presented are based on IMU2 processing, 
inertial trajectory estimates were obtained solving for state, attitude, and 
the 6 extended solution parameters previously described for all of the IMUs. 

Table m-1 shows the state vector solutions at the epoch time as well as an 
accuracy assessment. As can be seen, all 3 solutions compare favorably. 

The accuracy assessment was based on an ensemble of entry estimates and 
reflects a realistic judgment as to the accuracy with which the entry state is 
known. Formal statistics (1 0 ) as generated within ENTREE are generally 
several orders of magnitude smaller which is felt to be somewhat unrealistic. 

The state solutions obtained represent an "informat ion- only’’ solution that 
is, the results were completely determined from the tracking data content. 

The relatively large diagonal a priori covariance matrix used for the batch 
filter had virtually no effect on the solution. The data fits based on each of 
the three IMUs were essentially the same. The (RMSW) fits were 1. 14, 1. 15, 
and 1. 17 for IMUs 1, 2 and 3, respectively. This result shows that the data 
were fit to nearly lp in each case. This includes all the tracking data as well 
as the pseudo altimeter and pseudo Doppler data. 

Plots of selected planet relative and inertial parameters from the BET 
vs. time are shown in Figures m-la through IR-le, and vs. altitude in 
Figures IH-2a through IR-2e. These plots are based on the IMU2 estimate. 

The position and velocity are defined by: h, the geodetic altitude; 4^, the 
geodetic latitude; A , the longitude; V R , the planet relative velocity magnitude; 
y , the planet relative flight path angle; and ip R , the velocity vector heading 

relative to true North. Attitude angles, P R » $ R * and a R are the P lanet 
relative roll, sideslip, and angle of attack, respectively. The Euler angles, 
ip , d , and cp , are ordered yaw, pitch, and roll and define the attitude of the 
vehicle relative to a North- East-local vertical frame. The inertial velocity 
components relative to the same frame are given by u, v, and w, which are 
the North, East, and (positive) down components, respectively. Figures A-2 
and A-3 in Appendix A provide a graphical depiction of the attitude angles, 
position, and velocity components described above. 
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The estimate of the Shuttle position and velocily during runway rollout 
Is depicted in Figure in-3. Here the X-coordinate is measured along Runway 
23 from the surveyed runway threshold, positive in the direction of the Shuttle 
motion. Y is perpendicular to X in the horizontal plane, positive right as 
seen by the landing Shuttle. The altitude components are depicted in the 
bottom plots of Figure HI-3. Naturally, the actual terminal Shuttle velocities 
are zero post-stop, and the altitude of the c. g. above the runway during roll- 
out and under static conditions is approximately 4.8768 m (which is shown as 
a dashed line starting from nosewheel touchdown at t = 2317fo>. Also shown 
as dashed lines starting at t = 2368. 0 are the surveyed coordinate stop points 
(corrected for main wheel/center-of-gravity displacement) as measured 
following the flight: X = 4588 m; Y = -4. 4 m (F. O. E. D. Sketch No. 5120, 

Dryden Flight Research Center). 

The estimated stop position components are given in Table III-2. The 
estimated position at the stop time of 2368?0 was 15.2 in front of the surveyed 
stop point, 1. 2 m to the right, and 0.4 m high. The velocity difference esti- 
mates were all less than 0. 03 mps. The exceptional terminal altitude and 
velocity estimates are attributed to the processing of the pseudo altimeter and 
Doppler data (see Section II). The terminal state vector solutions for each of 
the 3 IMU-generated BETs are tabulated in Table III-2. 

Figures ni-4a through m-4j are the observation residual plots of all 
the measurement data processed in the generation of the BET associated 
with IMU2. Each page illustrates the data from a particular tracking station. 

The first plot shows the Guam S-band residuals. The next eight plots are the 
C-band residuals for PTPC, PPTC, VDBC, VDFC, VDSC, SNIC, FRCC, and 
EAFC, respectively. The radar types are noted thereon for each C-band station. 
The last figure contains residual plots for the three pseudo Doppler stations 
and altimeter observations. The left column on each figure shows the actual 
measurement residuals, O-C. The right column illustrates the weighted 
residuals, that is, the quotient of the actual residuals and the measurement 
weights. The computed means and standard deviations for each residual plot 
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are annotated thereon. Roundoff results in some of these quantities being 
displayed as absolute zeros. A weighted residual statistics summary is 
presented in Table III— 3. 

Generally speaking, the overall data fit is excellent. As can be seen 
from the residual plots, some slight signature trends remain, probably due 
to unmodeled error sources associated with the trackers and the IMUs. 
Nevertheless, with the exception of the range measurements from the PTPC 
station at the Point Pillar complex, all station residual statistics show means 
and standard deviations of less than , with most having a better than 1 Q fit. 

Table HI-3 also indicates that the residual spread and data fit are gen- 
erally independent of the dynamic data source. Most stations had either an 
all positive or all negative mean bias. Some were quite consistent in magni- 
tude. Note too that the pseudo altimeter had similar means and sigmas inde- 
pendent of the IMU used to generate the BET, whereas the pseudo Doppler data 
residual statistics for each IMU bore little resemblance to one another. 

Table HI-4 lists the IMU systematic error solutions associated with each 
of the inertial platforms. IMU1 yielded the smallest estimated accelerometer 
scale factor solutions. IMU3 yielded the smallest gyro drift bias estimates 
but the largest accelerometer scale factor error solutions. In general, the 
scale factor solutions showed the most consistency as the extended solve-for 
parameter set was varied. Indeed, the formal uncertainties associated with 
the scale factor solutions with all IMU modeled errors considered were gen- 
erally on the order of 50 - 100 ppm, indicative of a reasonably accurate 
estimate (the IMU specification accuracy as discussed in Appendix A is 100 ppm). 
On the other hand, the gyro drift bias solutions were very sensitive to nearly 
any change in the solution parameter set. Information only (i. e. , no a priori 
uncertainties) were 20 to 50 times larger than the gyro drift specification 
accuracies. There was insufficient information in the tracking data to obtain 
reliably accurate estimates of these parameters. 

Final atmosphere and atmosphere relative parameters are presented as 
Figs. ni-5a through III-5i. The atmosphere utilized was the Langley Atmos- 
pheric Information Retrieval file (LAIRS, USE8 dated October, 1981). 
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Figs. HI 5a through m-5d are plots of the temperature, pressure, density, 
and atmospheric wind profiles from this file. The winds are measured winds 
and are in general agreement with in situ determined winds as reported in 
Ref. 14. Also, additional measurements made at two California sites, 
Tehachapi and Wheeler Ridge, yielded similar wind profiles. The large planet 
relative side-slip angle excursions (~ 3 deg) shown in Figure RI-lc are due 
almost entirely to neglecting these winds in the attitude computation. 

Atmospheric relative velocity, flight path angle, and heading angle are 
shown in Fig. III-5e versus time. Air relative angle-of-attack and sideslip 
angle versus time are shown as Fig. ni-5f. Here it is shown that the air 
relative side-slip is within ±1.0 degree after inclusion of the atmospheric 
winds. This result is more reasonable and as anticipated based on STS-1 
measured spacecraft rudder deflections and lateral accelerations. Dynamic 
pressure and Mach No. time histories are shown as Figs. m-5g. Flight 
derived lift and drag coefficients as well as the L/D ratio are shown as 
Fig. III-5h. Also shown thereon are the flight derived side force coefficient 
versus time. Finally, flight derived pitching moment (C m ), yawing moment 
(C n ), and rolling moment (C^) coefficients are presented in Fig. m-5i. These 
air relative parameters are utilized by ACME investigators for post-flight 
assessments of the aerodynamic performance by comparing with preflight 
aerodynamic data base values. It is observed that the derived aerodynamic 
parameters do not stabilize until t ~ 700 sec due to the low signal to noise 
ratio of the measured rates and accelerations in the lowq environment. 
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Parameter 

Units 

IMU1 

IMU2 

IMU3 

1 a Accuracy 
Assessment 

V R 

km/sec 

7.41103 

7.41108 

7.41107 

l.E-4 

r R 

deg 

-1.1475205 

-1.1555853 

-1. 1530949 

4. E-3 

*R 

deg 

47.216922 

47.218146 

47.214843 

.01 

h 

km 

182.398 

182.994 

182.823 

0.250 


deg 

1.9323945 

1.9339547 

1.9333110 

1. E-3 

X 

deg 

140.76175 

140.76133 

140. 76203 

2. E-3 

o 

deg 

-7.4015553 

-7.4168490 

-7.3679519 


R 






8 

deg 

-1.4950769 

-1.5257547 

-1. 5227536 


R 






Or 

deg 

35. 548636 

35.592728 

35.585570 


rt 

* 

deg 

43.481720 

43. 494063 

43. 523341 

.08 

e 

deg 

34.255158 

34. 293573 

34.291767 

.02 

<p 

deg 

-8.9983262 

-9.0219117 

-8.9621916 

.05 

u 

km/sec 

5.0327 

5. 0327 

5.0330 


V 

km/sec 

5.4381 

5.4382 

5.4379 


w 

km/sec 

0.1484 

0. 1495 

0. 1491 



TABLE m-1 

STS-1 BET results at epoch using the tri-redundant IMUs 
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STATE VECTOR COMPONENT 
(RUNWAY COORDINATES) 

IMU1 

IMU2 

IMU3 

MEASURED END 
CONDITIONS 

x (km) 

4.6229 

4. 6037 

4.6000 

4.5884 

y (km) 

0. 0037 

-0. 0032 

-0.0064 

-0. 0044 

h (km) 

0.0051 

0. 0052 

0.0051 

I! 

x (mps) 

0.021 

0. 006 

0.021 

0.0 

y (mps) 

-0. 024 

-0. 018 

-0.021 

0.0 

h (mps) 

-0.018 

-0.027 

-0. 027 

0.0 


TABLE III-2 

BET terminal flight conditions from the tri-redundant IMUs for STS-1 
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II II II I lllllliMI! 


Station 


Weighted Mean, ju w 

Weighted Standard Deviation, ty w 

Data 

Type 

IMU1 

IMU2 

IMU3 

IMU1 

IMU2 

IMU3 

GWMS 

Range 

.06 

-.04 

.04 

.61 

.67 

.63 


Doppler 

.22 

-.33 

-.15 

1.00 

1.20 

1.00 


X- Angle 

-1.20 

1.21 

.50 

.70 

.70 

.57 


Y- Angle 

1.80 

.95 

.96 

.63 

.52 

.52 

PTPC 

Range 

-1.94 

-1.10 

-1.58 

.49 

.31 

.44 


Azimuth 

.62 

.61 

.68 

.41 

.41 

.41 


Elevation 

.80 

.94 

.82 

.28 

.30 

.29 

VDBC 

Range 

-1. 01 

-.67 

-.54 

.67 

..94 

.96 


Azimuth 

-.22 

-.25 

-. 33 

.47 

.47 

. 54 


Elevation 

.26 

.44 

.32 

.31 

.31 

.32 

VDFC 

Range 

-1.37 

-1.01 

-.89 

1. 12 

1.50 

1.51 


Azimuth 

-.06 

-.09 

-.17 

.86 

.88 

• 95 


Elevation 

.23 

.37 

.25 

.71 

.78 

.76 

VDSC 

Range 

-.16 

.17 

.28 

.65 

1.00 

1.03 


Azimuth 

-.23 

-.26 

-.34 

.85 

.86 

. 95 


Elevation 

-. 14 

0.0 

-.13 

,74 

.82 

.77 

SNIC 

Range 

.57 

.53 

.50 

.90 

.99 

.97 


Azimuth 

-1.69 

-1.68 

-1.66 

.73 

.72 

. 74 


Elevation 

-. 07 

-.01 

-.03 

.89 

.90 

.89 

FRCC 

Range 

-.63 

-. 90 

-1. 04 

1.17 

1.13 

1.12 


Azimuth 

.15 

.32 

,27 

.94 

1. 15 

1. 16 


Elevation 

.76 

.86 

.80 

.85 

.86 

.80 

EAFC 

Range 

-.04 

-.24 

-.36 

1.26 

1.11 

1.12 


Azimuth 

-. 05 

.11 

.08 

1.08 

1.22 

1.23 


Elevation 

.62 

.75 

.77 

1.04 

1.05 

1. 19 

PPTC 

Range 

-3. 05 

-2.16 

-2.51 

.63 

.51 

.56 


Azimuth 

.04 

.07 

.27 

.27 

.32 

.38 


Elevation 

-.70 

-.52 

-.79 

.53 

.53 

. 46 

Pseudo 

Altimeter 

-.37 

-.50 

-.38 

.14 

.21 

.21 

Pseudo 

Doppler# 1 

.03 

-.07 

.18 

.58 

1.25 

.74 


Doppler#2 

-.03 

.56 

.83 

.95 

1.09 

1. 31 


Doppler#3 

1.01 

2.15 

1.58 

.62 

.90 

. 61 


TABLE m-3 

Weighted residual statistics summary for STS-1 
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IMU1 

IMU2 

IMU3 

X-gyro drift 

-0. 146 deg/hr 

-0.092 deg/hr 

+0, 050 deg/hr 

Y-gyro drift 

“0. 051 deg/hr 

+0 , 110 deg/hr 

-0. 021 deg/hr 

Z-gyro drift 

-0. 012 deg/hr 

+0, 096 deg/hr 

+0. 020 deg/hr 

X-accelerometer 
scale factor 

-8 ppm 

56 ppm 

| 

193 ppm 

Y-accelerometer 
scale factor 

-16 ppm 

190 ppm 

162 ppm 

Z -a c celer ometer 
scale factor 

13 ppm 

-64 ppm 

-144 ppm 


TABLE HI-4 

IMU parameter estimates for STS-1 
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Figure m-2e. 
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Figure m-4a. STS-1 final Guam S-band 
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Figure m-4b. STS-1 final Pt. Pillar (PTPC/FPQ-6) residuals versus time 
from epoch 
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Figure Itl-4c. STS-1 final Pt. Pillar (PPTC/FPS-16) residuals versus time 
from epoch 
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Figure m- 4 d. STS -1 final Vandenberg (VDBC/TPQ- 18 ) residuals versus 
time from epoch 
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Figure m-4e. STS-1 final Vandenberg (VDFC/FPS-16) residuals versus 
time from epoch 
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Figure m-4f. STS-1 final Vandenberg (VDSC/FPS-16) residuals versus 
time from epoch 
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Figure EI-4g. STS-1 final St. Nicolas Island (SNIC/FPS-16) residuals versus 
time from epoch 
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Figure IH-4h. STS-1 final NASA Dryden (FRCC/FPS-16) residuals versus 
time from epoch 
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Figure m-4i. STS-1 final Edwards (EAFC/FPS-16) residuals versus 
time from epoch 
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Figure HI-5f. STS-1 BET atmospheric relative angle-of-attack and 
side-slip angle versus time from epoch 
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IV. Summary 

The STS-1 Space Shuttle re-entry trajectory has been successfully re- 
constructed using a weighted least squares batch filter algorithm. Dynamic 
data derived from the onboard Inertial Measurement Units (IMU) were used 
to propagate the state vector. Tracking data from eight California based 
C-band radar stations and the S-band tracking station at Guam were pro- 
cessed in the BET generation. The Guam data in particular were instru- 
mental in anchoring the position and velocity estimates at ~ 183 km altitude. 
Likewise, the pseudo altimeter and pseudo Doppler data processed during and 
post rollout significantly improved the estimation accuracy during the terminal 
portion of the trajectory. 

Examination of the BET output demonstrated that the STS-1 re-entry 
trajectory was quite similar to the pre-mission nominal flight profile. IMU 
to IMU comparisons, and IMU systematic error solutions indicated nominal 
platform performance. Processing selected data from all available tracking 
stations resulted in an approximate 1 a overall RMSW fit for each of the 
3 IMU determined BETs, thus generating confidence in the accuracy of the 
estimation. In summary, the important in-plane entry parameters ( V, y , h) 
were determinable (1 c) to 0.01 mps, 0.004 deg, and 250 m, respectively. 
Spacecraft attitude accuracies at epoch of 0. 08 deg, 0. 02 deg, and 0. 05 deg 
are estimated for the inertial Euler angles $ , 6 , and <p , respectively. 
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APPENDIX A 

Discussion of the BET Generation Process 
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This Appendix is presented to provide for a general discussion of the 
data pre-processing required to enable the generation of a BET. Tracking 
data and dynamic data pre-processing requirements are addressed. A soft- 
ware overview is shown as Figure A-l. Table A-l presents a list of acro- 
nyms for the software referred to herein. The overall ENTREE software 
Bystem is summarized to show the data flow between receipt of data to gen- 
eration of the final BET for the user community. Shuttle specific pre- 
processing requirements developed by AMA, Inc. under the subject contract 
to satisfy the ENTREE software are addressed. Pre-processing peculiar 
to the STS-1 flight are addressed in the text of the report. The output pro- 
duct from ENTREE is an inertial BET. The final product, as shown in 
Figure A-l, combines the ENTREE output with the best available atmos- 
phere information (including winds). The atmosphere is provided by LaRC, 
with contractual help from the Space Systems Division of Computer Sciences 
Corporation, in the form of a Langley Atmospheric Information Retrieval 
System file. This atmosphere is developed from a combination of measure- 
ments and models as discussed in Ref. 5 and is translated in time and space 
to conform to the ground track and vertical profile of the BET. These data 
permit the computation of the required air relative parameters and, along 
with the measured accelerations, rates, and Shuttle mass properties, enables 
computation of flight derived aerodynamic force and moment coefficients. 

A. 1 ENTREE Software Description 

The major estimation software, ENTREE (Ref. 8), was initially 
developed by the Computer Sciences Corporation under Contract NAS1-15663 
for LaRC. AMA, under the subject contract, has had considerable involve- 
ment in checkout and modifications/additions to this software. The software 
requires body -fixed (strapped-down) dynamic measurements for use in the 
six-degrees-of-freedom equations of motion for spacecraft prediction. Body 
axes conventions for the angular rates and linear accelerations conform to the 
usual aerodynamicists' definitions as depicted in Fig. A-2. A fourth order 
fixed step size Runge-Kutta integration algorithm is utilized. Definition of 
the variables utilized in the software can best be described by referring to 
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Figures A-3 a,b. Figure A-3a shows the planet model, position, and velocity 
parameters. The altitude corresponds to an altitude above an oblate spheroid 
which conforms to the Fischer model. Longitude, X , is defined as positive 
Eastward from Greenwich. Inertial velocity components, u, v, and w, are 
geocentrically oriented to local North, East, and vertical (downward). The 
velocity heading angle, flr , is defined positive clockwise from North and the 
flight path angle, y , is defined positive above the geocentric horizon. Space- 
craft attitude parameters are shown as Figure A-3b. The velocity relative 
parameters are: 0 , roll with respect to the velocity vector (positive right 
wing down); 0 , side-slip angle (positive nose left); and a, the angle-of-attack 
positive (nose up). Geocentrically oriented Euler angles are also utilized. The se 

quence Is yaw, $ , pitch , 0 , and roll , <p , and orients the vehicle body axes 
to the local vertical system. Though not shown in the schematic, a software 
utility, TRANS, has been developed to compute the required ENTREE state 
variables from the initial state estimate in the inertial 1950. 0 Mean Equator 
and Equinox (M50) system. Also, based on this M50 state and interpolated 
IMU measurements at epoch, initial attitude estimates are generated therein. 

Batch weighted least squares and sequential Kalman filtering algorithms 
can be selected on option for the estimator. A weighted least squares batch 
filter is employed to obtain the best estimate based on the observations pro- 
cessed. 

Potential observables which can be selected on option (see Refs. 8, 13, and 15) 
are: 

C-band Range, Azimuth, and Elevation 
S-band Range, Doppler, X-angle, and Y-angle 
Tacan Range, and Bearing angle 
Altimeter 

Microwave Scanning Beam Range, Azimuth, and Wedge angle. 

Of particular importance for Shuttle are the C-and S-band observables. 

Tacan accuracy, relative to these radars, and MSBLS timing staleness in 
the down-list do not warrant use of these observables. 


60 



A. 2 Tracking data pre-processing 

Two software utilities have been developed, PREOBS and 
OBEDIT , to employ the external observations in ENTREE (see Figure A-l). 
PREOBS reads the tracking data files from several sources, i.e. , GSFC, 

JSC, and recorded 01 data. These data are transmitted to LaRC and con- 
verted by the Orbit er Experiments (OEX) Data Manager to be compatible with 
the LaRC computer system. 

The GSFC input as shown represents the primary source for high speed 
S-band tracking prior to the entry interface. These GSFC data were obtained 
through special arrangements with LaRC. These data are playback data. The 
necessity for the high rate data is as follows. The ENTREE program uses 
a modified formulation of an instantaneous range rate computation for Doppler 
frequency shift. Since the S-band Doppler measurement is accumulated cycles 
over a time interval (count time) and must be converted to frequency, an in- 
stantaneous formulation requires a very small count time for accuracy. Prior 
to entry interface the real time data are transmitted to the JSC at a 10 second 
rate which is unacceptably large in terms of count time. 

Range, Doppler, X-angle, and Y-angle measurements are all included on 
the GSFC file. Low rate S-band data are also contained on the JSC tracking 
file prior to the entry interface. Use is made of these data to check on time 
tags for the high rate (playback) data from GSFC. The principal measure- 
ments taken from the JSC tracking data file are the C-band tracking data be- 
tween end of communications blackout and touchdown. The C-band measure- 
ments (Range, Azimuth, Elevation) provided on the JSC file are in units 
compatible with ENTREE and require no units conversions or calibrations. 
S-band X and Y-angle measurements obtained from the JSC file are in units 

compatible with ENTREE. Those obtained from the GSFC file are converted 

from angle units (where one unit is a specified number of degrees) to radians. 

S-band ranging measurements are in fact round trip light time measure- 
ments. As such they must be calibrated for timing delays occurring at both 
the station and the spacecraft. For Shuttle, S-band ranging measurements are 
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calibrated "on site" for station delays but not the spacecraft delay. The 
signal turn around delay in the spacecraft S-band ranging transponder 
varies slightly over a station pass. This transponder delay is assumed 
constant, however, and is subtracted from each S-band ranging measure- 
ment. The value of the transponder delay is provided by the JSC. The S-band 
ranging measurements on the GSFC file are in units of round trip light time 
and are converted to average slant range. The S-band ranging measurements 
on the JSC file have already been converted to average slant range. In either 
case, the ranging measurement is "calibrated" by decreasing its value by the 
range equivalent of the transponder delay* 

S-band Doppler data from either GSFC or JSC are provided as counted 
cycles. Doppler frequency is obtained by differencing the counter 
readings, dividing by the count time and then subtracting the frequency bias. 

The resulting "measurement", which may be thought of as average slant range 
rate over the count interval, is time-tagged at the midpoint of the count interval 
to better approximate instantaneous slant range rat*;. 

On option, the alternate datatypes, TACAN, MSBLS, and altimeter, are 
obtained from the spacecraft recorded data as separate files. At present, no 
use is made of these data for entry reconstruction though pseudo altimeter 
measurements were processed to improve the BET during rollout for STS-1. 

Software PREOBS reads the tracking data files and merges and orders 
by time and station all the data types for ENTREE processing. During the 
estimation process blunder points can be rejected within ENTREE, either by 
sigma rejection or elevation masking. Another tracking data processor, 

OB EDIT, maybe used as a preprocessor but it is really an "in-line" processor. 
OB EDIT is used for time deletion of selected measurements on the ENTREE 
input tracking data file. The "selected" measurements are either isolated 
blunder points or a group of measurements over a time interval. An examin- 
ation of post-fit residuals is used in determining which data are to be deleted 
from the tracking file prior to the next ENTREE estimation run. 
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A. 3 Dynamic Data 


There are four potential sources of dynamic data available 
for use in ENTREE. There are the strapped-down measurements from 
the Aerodynamic Coefficient Identification Package (ACIP)^ and the mea- 
surements from the tri -redundant IMUs. Though the ACIP measurements 
satisfy the ENTREE strapped-down requirements, pre-flight test results 

(Ref. 16) indicated that these data were not of sufficient accuracy to utilize 
in the BET generation. (The ACIP data are of sufficient accuracy to 
extract aerodynamic coefficients and, because of the high frequency (~ 170 Hz) 
of the measurements, are utilized by MMLE investigators to extract stability 
derivatives and aerodynamic control surface effectiveness). Therefore, this 
discussion focuses on the utilization of the tri-redundant IMUs to satisfy the 
ENTREE interface. 

EMU data are obtained via the JSC. These data are also converted by the 
OEX Data Manager for LaRC use. IMU pre-processing requirements are 
two-fold. First, due to the redundant nature of the IMUs, comparisons must 
be made to define, at least on a relative basis, the performance of the tri- 
redundant set. Secondly, pre-processing to emulate the required strapped- 
down measurements is required. 

The tri-redundant IMUs are gimballed inertial platforms whose orien- 
tations are skewed with respect to one another and are located at the navigation 

. ; - ( 2 ) 

base in the nose of the Shuttle vehicle. The 1 <7 accuracy specifications for 
these units are defined in Ref. 17 and listed here: 

accelerometer bias: 50 jU g (10 ^g) 
accelerometer scale factor : 100 ppm 
gyro drift bias; .035 deg/hr (.022 deg/hr) 
gyro g-sensitive drift bias: . 025 deg/hr/g 
initial platform misalignments: (80 sec) 


1 The simplified schematic. Figure A-l, does not show any pre-processing 
refinements to utilize the ACIP data in ENTREE. It should be understood 
that, at a minimum, comparisons of ACIP measurements with derived IMU 
body axis data are required. 

^Numbers in parentheses presume pre-deorbit calibrations and star tracker 
alignment. 



Additionally, the IMU accumulated velocity output as measured by the accel- 
erometers is quantized to 1 cm/sec. Likewise, the gyro gimbal resolver 
output, the ultimate source of the platform to outer roll quaternion, is 
quantized to multiples of 20 sec. 

The output of each IMU consists of the 3 components of accumulated 
sensed velocity, expressed in M50 coordinates, and the 4 components of the 
platform to outer roll quaternion. This output is available from the real time 
telemetry data and is simultaneously recorded onboard. Because the IMU 
output data rate differs from the downlist (D/L) sequencer data rate, the most 
frequent IMU output (6. 25 Hz) is not time tagged and use of these data 
was not considered. However, time tags associated with the velocity (and 
quaternion) components are stored and recorded within the D/L frame at 
approximately 1 Hz in order to insure data homogeneity. These data are not 
at a uniform rate. For example, the 4 quaternion components of all 3 IMUs 
are simultaneously output at a 0. 96 second rate. With a 1. 0 second D/L rate, 
each quaternion output record on the T/M tape differs In time from the previous 
record by 0. 96 seconds, except for every 24th record which jumps to 1. 92 sec 
when two quaternion output records fall within the same D/L frame and the first 
is overwritten. The same holds true for the velocity components of the IMUs 
(although time tagged different from the quaternion data) with the exception of an 
output rate change from 0. 96 seconds to 0. 16 seconds starting at the initiali- 
zation of the entry guidance mode 5 minutes prior to entry interface. This 
change results in an input velocity record spacing of 0. 96, 0. 96, 0. 96, 1. 12, 
0.96, 0.96, 0.96, 1.12 (seconds), etc., thereafter. 

Selection of the best IMU for use in ENTREE is of utmost importance. 

A procedure has been established to compare independently the gyro and 
accelerometer performance of each IMU versus the remaining two as well as 
combinations of the measurements from the various sets. This procedure, 
and STS-1 results, are discussed in Ref. 11 and briefly summarized here. 

Figure A-l shows the software flow to enable the mutual comparisons, speci- 
fically the utilities PREVEL, ABSATT and CALIBRT. PREVEL provides 
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a measure of accelerometer performance by comparing M50 velocity measurements. 
These comparisons are not independent of gyro performance since the orienta- 
tion of each platform with respect to the inertfal frame is assumed absolutely 
known. ABSATT provides for a measure of gyro performance by comparing 
inertially referenced Euler angles as suggested independently by the tri- 
redundant set. Finally, the software utility, CALIBRT, determines first order 
calibrations, e.g. , accelerometer scale factors, gyro drifts, accelerometer 
biases, of each IMU with respect to some selected fiducial reference set. 

The major software required to satisfy the ENTREE interface is 
PREIMU. PREIMU, operating from the reformatted, edited, file generated 
by PRETM, derives the equivalent spacecraft rates and accelerations in the 
platform axes. Transformation to body axes and accommodation of sensor 
locations with respect to the Shuttle center-of-gravity are done internal to 
ENTREE. PREIMU processing of the IMU data into a form compatible for 
dynamic data input to ENTREE is described in detail in Reference 18. In 
summary, the M50 velocities are spline fitted and differentiated to yield an 
acceleration time history (which, when integrated, yields the original velocity 
history by definition) at a user defined rate with any data gaps filled, if re- 
quired. The accelerations are rotated to platform coordinates using the 
REFSMMATs (see Table B-2 in Appendix B) and stored on the ENTREE input 
dynamic data file. The platform to outer roll quaternion information is com- 
bined with pad loaded navigation base to body and navigation base to outer roll 
transformation matrices to produce a set of platform to body Euler angles (or 
quaternions). These angles (quaternions) can then be spline fitted and differen- 
tiated to yield Euler angle rates (quaternion rates) at the same times as the 
acceleration data. The transformation to angular rates about the IMU X, Y, 
and Z axes is then straightforward. These rates are also stored on the ENTREE 
input dynamic data file, along with the platform to body Euler angles (or quater- 
nions). These 11 element data records (time, platform attitude rates (3), plat- 
form accelerations (3), and quaternions (4) (or Euler angles (3) plus a flag (1) ) 
provide the necessary information for ENTREE to solve for systematic IMU 
errors in the platform coordinate system as well as integrate the equations of 
motion in the strapped-down coordinate system. 
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As just described, the preprocessor program has the option of appending 
the platform to body attitude information to the dynamic data input file in the 
form of either quaternions or Euler angles. Furthermore, two of the 12 
potential Euler angle sequences are programmed as options, with the before- 
mentioned flag value signifying the sequence chosen. Each option has potential 
disadvantages. The differentiated quaternion data cannot be guaranteed to 
yield orthonormal transformations, while an Euler angle sequence could con- 
ceivably result in a singularity condition at a certain platform to body attitude. 
As it turned out, the Euler angle sequence chosen for the STS-1 post flight 
processing did not encounter any singularities. 

As stated previously, the manipulations required to pre-process the 
IMU data result from the use of an inertial instrument’s data in a strap-down 
formulation. The use of the Aerodynamic Coefficient Identification Package 
(ACIP) with its body mounted linear accelerometers and rate gyros would be 
a natural for input data. Unfortunately, the accuracy specifications associated 
with the ACIP preclude its use for BET generation. 
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acronym 

ABSATT 

CALIBRT 

ENTREE 
MMLE 
NEWBET 
OB EDIT 
PREIMU 

PREOBS 

PRETM 

PREVEL 

TRANS 


FUNCTION 

Absolute IMU attitude measurement comparison 
software 

IMU calibration software for first order perfor- 
mance comparisons 

Entry Trajectory Reconstruction Software 

Modified Maximum Likelihood Estimator 

Software to merge inertial BET and atmosphere 

Observation data editor 

Cubic spline processor to derive spacecraft rates 
and accelerations from IMU measurements 

Software to pre-process observation data from 
available sources 

Software to pre-process and edit IMU data 

IMU accelerometer performance comparison 
software for M50 velocity measurements 

Software to transform inertial M50 initial state 
estimates to ENTREE coordinates 


TABLE A-l 
Software Acronyms 
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Figure A-l. Schematic of software/data interfaces required to 
generate BET 
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APPENDIX B 

STS-1 MISSION SPECIFIC INPUT DATA 
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This Appendix contains STS-1 mission specific input data required to 
generate the BET. Table B-l presents the station characteristics which 
includes type, internal numbering system and associated acronym utilized, 
the best location set for metric data processing, station frequency and radar 
mount if applicable, index of refraction based on the mean monthly average 
for April, 1981, and the atmospheric scale height utilized in the refraction 
modelling. Table B-2 presents the relevant attitude matrices required to 
process the IMU measurements to derive body axis data. Table B-3 lists 
the elements of the a priori diagonal covariance matrix used in the batch 
solution. Finally, Table B-4 presents the inputs utilized for the planet model, 
runway location, IMU location with respect to the Shuttle center-of-gravity, 
and mass properties and associated aerodynamic reference parameters re- 
quired to compute the in-flight aerodynamic force and moment coefficients. 
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Station locations and refraction data for STS-1 data processing 




"REFSMMAT" MATRICES (M50 TO 

PLATFORM! 


(IMU1) 


-0.79266172 

-0.57519790 

-0. 20207602 

-0. 44474863 

+0. 77226827 

-0.45365167 

+0.41699673 

-0* 26971874 

-0. 86796603 


(TMU2) 


+0. 39075335 

-0.88816900 

-0.24179873 

+0. 73866644 

+0. 45929717 

-0.49337438 

+0. 54925762 

+0. 01417890 

+0. 83553258 


(IMU3) 


-0. 11996126 

-0. 55783236 

+0. 82123822 

+0. 22795183 

+0.78963381 

+0. 56966257 

-0. 96625435 

+0.25553983 

+0. 03243346 

NAV BASE TO OUTER ROLL Transformation Matrices 


(IMU1) 


0. 99999938 

-3. 9528892E-4 

-1. 0380259E-3 

3. 9528871E-4 

0. 99999992 

-4. 1032019E-7 

1. 0380260E-3 

0.0 

0. 99999946 


(IMU2) 


0. 9999979 

1. 9876E-3 

-5. 693E-4 

-1. 9876E-3 

0.9999980 

-6. 3E-6 

5.693E-4 

7. 5E-6 

0. 9999998 


(IMU3) 


0. 9999934 

3. 7771341E-3 

-3. 9463471 E-4 

-3. 7771382E-3 

0. 99999338 

-1. 0690797E-5 

3. 9459304E-4 

1. 1723317E-5 

1.0 

NAV BASE TO BODY (all IMUs! 

0.9829565 

4. 363323E-4 

-0. 1838379 

-4. 529508E-4 

0. 9999999 

-4. 84048 E-5 

0. 1838379 

1. 308493E-4 

0. 9829566 


TABLE B-2 


STS-1 Attitude transformation matrices required for IMU processing 




<p 

platform drift (each axis) = 0. 083 deg/hr 
accelerometer scale factor (each axis) = 400 ppm 




Planet Parameters 


Physical Model 


Polar Radius: 
Equatorial Radius: 
Rotational Rate: 


Gravity Model 


Central term, 


12 * 


J 4 : 

c 22 : 

S 22 : 


6356. 784284 km 
6378. 166 km 
, 7292115147E-4 rad/sec 


, 398601999995E15 m 3 /sec 2 


. 10827E-2 
. 256E-5 
. 158E-5 
. 157E-5 
. 897E-6 


Runway 23 Location: 


Altitude: 

Geodetic Latitude: 

Longitude: 

Azimuth: 


635. 8128 m (above ellipsoid) 
34. 966397 deg 
242. 180352E deg 
244.413472 deg 


Location of IMU relative to center-of-gravitv in Body coordinates 
(Assumed constant during Entry) 


X R 

Zt> 


17. 0688 m 
0. 0 m 
-1.2192 m 


STS-1 mass properties and aerodynamic reference parameters 


Weight 

Reference Area 

Span 

Chord 


Moments and products of inertia: 


89930.448 kg 
249. 909 m 2 
23. 792 m 
12. 060 m 


I ^ 1213866 kg-m 2 

Iyy 9378654 kg-m 2 

I zz 9759518 kg-m 2 

\ z 228209 kg-m 2 

Ixy 6136 kg— m 2 

Iy Z 2972 kg-m 2 
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TABLE B-4 

Planet and spacecraft data used for STS-1 BET generation 



APPENDIX C 


LISTING OF STS-1 BET PARAMETERS 



This Appendix is presented to provide a listing of the actual BET 
parameters at a reasonable spacing. The listing was generated from a 
permanent file (METBET1 under user catalog, UN = 274885C) which is 
the metric equivalent to STS1BET, that version in English units widely 
used by the user community at LaRC and the various other NASA agencies, 
including the AFFTC at Edwards and Rockwell personnel. Alphanumeric 
definition of the variables and units utilized are as defined in Ref. 7 and as 
noted on the listing of the header record. Above ~ 30 km, the data are 
presented at 50 sec intervals. The remainder of the data are given at a 

5 sec spacing. Both files, METBET1 and STS1BET, are actually written 
at 1 sec spacing. 
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